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Abstract 

The main objective of this work is to explore the evolution in the structure of 
the quark-anti-quark bound states in going down in the chirally restored phase 
from the so-called "zero binding points" T^b to the full (unquenched) QCD critical 
temperature Tc at which the Nambu-Goldstone and Wigner-Weyl modes meet. In 
doing this, we adopt the idea recently introduced by Shuryak and Zahed for charmed 
cc, light-quark qq mesons vr, cr, p, Ai and gluons that at Tzb, the quark-anti-quark 
scattering length goes through oo at which conformal invariance is restored, thereby 
transforming the matter into a near perfect fluid behaving hydro dynamically, as 
found at RHIC. We show that the binding of these states is accomplished by the 
combination of (i) the color Coulomb interaction, (ii) the relativistic effects, and 
(iii) the interaction induced by the instanton-anti-instanton molecules. The spin-spin 
forces turned out to be small. While near Tzb all mesons are large-size nonrelativistic 
objects bound by Coulomb attraction, near Tc they get much more tightly bound, 
with many-body collective interactions becoming important and making the a and 
vr masses approach zero (in the chiral limit). The wave function at the origin grows 
strongly with binding, and the near-local four-Fermi interactions induced by the 
instanton molecules play an increasingly more important role as the temperature 
moves downward toward Tc- 
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1 Introduction 



1.1 qq bound states above 

The concept that hadronic states may survive in the high temperature phase 
of QCD, the quark-gluon plasma, has been known for some time. In particular, 
it was explored by Brown et al.[l,2]. The properties of (degenerate) tt and a 
resonances above Tc in the context of the NJL model was discussed earlier 
by Hatsuda and Kunihiro[3], and in the instanton liquid model by Schafer 
and Shuryak [4]. Recently, lattice calculations [5,6] have shown that, contrary 
to the original suggestion by Matsui and Satz [7], the lowest charmonium 
states J/V', remain bound well above T^.. The estimates of the zero binding 
temperature for charmonium Tj/^ is now limited to the interval 2Tc > Tj/^ > 
1.6Tc, where ~ 270 MeV is that for quenched QCD. Similar results for light 
quark mesons exist but are less quantitative at the moment. However since 
the "quasiparticle" masses close to are large, they must be similar to those 
for charmonium states. 

In the chiral limit ^ all states above the chiral restoration go into chiral 
multiplets. For quark quasiparticles this is also true, but although the chirality 
is conserved during their propagation, they are not massless and move slowly 
near where their "chiral mass" m — E{p — > 0) is large (~ 1 GeV). 

RHIC experiments have found that hot/dense matter at temperatures above 
the critical value Tc{unquenched) ~ HQ MeV is not a weakly interacting gas 
of quasiparticles, as was widely expected. Instead, RHIC data have demon- 
strated the existence of very robust collective flow phenomena, well described 
by ideal hydrodynamics. Most decisive in reaching this conclusion was the 
early measurement of the elliptic flow which showed that equilibration in the 
new state of matter above Tc set in in a time < 1 fm/c [8]. Furthermore, the 
first viscosity estimates [9] show surprisingly low values, suggesting that this 
matter is the most perfect liquid known. Indeed, the ratio of shear viscosity 
coefficient to the entropy is only t]/ s ~ 0.1, two orders of magnitude less than 
for water. Furthermore, it is comparable to predictions in the infinite coupling 
limit (for N=A SUSY YM theory) r]/s = I/Att [10], perhaps the lowest value 
possible. 

Shuryak and Zahed[ll] (hereafter referred to as SZ whenever unambiguous) 
have recently connected these two issues together. They have suggested that 
large rescattering cross sections apparently present in hot matter at RHIC are 

^ In most of what follows, we will be ignoring the effect of light-quark masses, 
unless mentioned otherwise. The "quasiparticle mass" that we shall refer to in what 
follows is a chirally invariant object called "chiral mass." 
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generated by resonances near the zero-binding lines. Indeed, at the point of 
zero binding the scattering length a of the two constituents goes to cxd and 
this provides low viscosity. This phenomenon is analogous to the elliptic flow 
observed in the expansion of trapped ^Li atoms rendered possible by tuning 
the scattering length to very large values via a Feshbach resonance [12]. 

Near the zero-binding points, to be denoted by T^t, introduced by SZ the 
binding is small and thus the description of the system can be simple and 
nonrelativistic. The binding comes about chiefly from the attractive Coulomb 
color electric field, as evidenced in lattice gauge calculation of Karsch and 
collaborators [5, 13], and Asakawa and Hatsuda[6], as we shall detail. The in- 
stanton molecule interactions, which we describe below, are less important at 
these high temperatures (T ~ 400 MeV). All this changes as one attempts 
(as we show below) to discuss the more deeply bound states just above Tc 
(unquenched). 

In another work [14], Shuryak and Zahed have also found sets of highly rela- 
tivistic bound light states in the strongly coupled iV=4 supersymmetric Yang- 
Mills theory at finite temperature (already mentioned above in respect to vis- 
cosity). They suggested that the very strong Coulomb attraction can be bal- 
anced by high angular momentum, producing light states with masses m ~ T. 
Furthermore, the density of such states remains constant at arbitrarily large 
coupUng. They argued that in this theory a transition from weak to strong 
coupling basically implies a smooth transition from a gas of quasiparticles to a 
gas of "dimers" ^ , without a phase transition. This is an important part of the 
overall emerging picture, relating strong coupling, viscosity and light bound 
states. 

In this work we wish to construct the link between the chirally broken state 
of hadronic matter below Tc (unquenched) and the chirally restored mesonic, 
glueball state above Tc- Our objective is to understand and to work out in 
detail what exactly happens with hadronic states at temperatures between 
Tc and T^h- One important new point is that these chirally restored hadrons 
are so small that the color charges are locked into the hadrons at such short 
distances (< 0.5 fm) that the Debye screening is unimportant. This is strictly 
true at T ^ T^, where there is very little free charge. In this temperature range 
the nonrelativistic treatment of SZ should be changed to a relativistic one. 

The relativistic current- current interaction, ultimately related with the clas- 
sical Ampere law, is about as important as the Coulomb one, effectively dou- 
bling the attraction (see section 2.1). We also found that the spin-spin forces 
discussed in 2.2 are truly negligible. In effect, with the help of the instanton 

^ In the large number of color limit considered, those are dominated by bound 
colored states, not colorless mesons that are important in the regime we are consid- 
ering. 
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molecule interaction, one can get the bound quark-anti-quark states down in 
energy, reaching the massless a and tt at T^., so that a smooth transition can 
be made with the chiral breaking at T < T^. 

The nonperturbative interaction from the instanton molecules becomes very 
important ^ . Let us remind the reader of the history of the issue. The nonper- 
turbative gluon condensate, contributing to the dilatational charge or trace 
of the stress tensor T^^ = e — Sp, is not melted at Tc (unquenched). In fact 
more than half of the vacuum gluon condensate value remains at T right above 
Tc- the hard glue or epoxy which explicitly breaks scale invariance but is un- 
connected with hadronic masses. The rate at which the epoxy is melted can 
be measured by lattice gauge simulations, and this tells us the rate at which 
the instanton molecules are broken up with increasing temperature. We will 
discuss this further in section 1.3. 

As argued by Ilgcnfritz and Shuryak [16] (for further references see the review 
[15]), this phenomenon can be explained by breaking of the instanton ensemble 
into instanton molecules with zero topological charge. Such molecules generate 
a new form of effective multi-fermion effective interactions similar to the orig- 
nal NJL model, see details in section 2.3. Brown ct al.[17] (denoted as BGLR 
below) obtained the interaction induced by the instanton molecules above 
by continuing upwards the Nambu-Jona-Lasinio description from below Tc- 

Our present discussion of mesonic bound states should not be confused with 
quasi-hadronic states found in early lattice calculations[18] for quarks and 
anti-quarks propagating in the space-like direction. Their spectrum, known 
as "screening masses" is generated mostly by "dynamical confinement" of 
the spatial Wilson loop which is a nonperturbative phenomenon seen via the 
lattice calculations. Similar effects will be given here by the instanton molecule 
interaction. We will briefly discuss it in section 1.4. 

1.2 Quasiparticles and their masses atT > Tc 

It is well known that the expectation value of the Polyakov line | (L) | goes to 
zero at Tc, indicating an infinite quark mass below T^; i.e., confinement. In the 
deconfined phase we discuss this mass is finite, and as T grows it is expected 
to decrease to some minimal value, before growing perturbatively as M ~ gT 
at large T. 

Chirality of the quarks is a good quantum number above T^. Now the chirally 
restored wave functions have good helicity, -|- or — . Also it can be seen from 

^ Although we do not discuss T < in this work, we still mention that in this 
region the instanton-induced effects seem to become dominant, see [15] for review. 
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eq. (9) that each wave function has good chirality. a ■ p — ±p ip. The 
hehcities and chirahties use up the 4 components in the Dirac wave function. 
The fermion modes with equal (opposite) hehcities and chirahties are called 
"quark quasiparticles" ("plasmino" modes). 

There are detailed results on their dispersion relations from pQCD. Specifically 
for quark quasiparticles, Weldon[19] obtained the dispersion relation 



which can be reasonably well approximated (to within 10%) by pi ^ +p^ 
where = The perturbative gluon mass is = g'^T'^{Nj2 + 

Nf/6)/2 [20]. 

If the exact dispersion relations for quasiparticles are known, the 2-body inter- 
action can be introduced by standard substitution of frequency and momenta 
by the covariant derivatives with the potentials. 

Unfortunately, the actual lattice data about the masses are very fragmented, 
and are (to our knowledge) only available from quenched calculations. The 
Polyakov loop expectation value \ {L)\{T) is plotted in Asakawa et al.: at each 
Nt it starts to deviate from zero at T^. ~ 270 MeV, the critical temperature for 
quenched calculations. However, one has to renormalize the mass, removing 
the linear divergent part of the point charge mass ~ 1/a before extracting any 
numbers, which needs more data. 

There are some results on quasiparticle dispersion relations from Petreczky 
et al. [13], who find in (quenched) calculation that they are all consistent 
with the usual dispersion relation ut'^ = k"^ + w? with the masses rUq/T = 
3.9 ± 0.2, rrig/T = 3.4 ± 0.3 or rUq ~ 1.6 GeV , rrig - 1.4 GeV for T = 
1.5Tc = 405 MeV where we take Tc = 270 MeV. These quarks (and gluons) 
are about as massive as charmed quarks. Although obviously these masses 
will be reduced in unquenched calculations, probably roughly by the ratio of 
critical temperatures or so, it seems inevitable that at T ^ Tc the quarks are 
so massive ^ that their Boltzmann factors are small, inadequate to furnish 
substantial presence of free quarks and gluons. The light bound state would 
then be primary sources of the pressure, as at low T. 

Petreczky et al.[13] omitted the first three values (at lowest temperatures) 
from their fit of uJq{p). Their curve of uj'^{p)/T'^ intercepts the ordinate at a 



^ We will argue for a somewhat lower values below, which would perhaps arise 
from the unquenched calculations. Let us also remember that those are the "chiral 
masses" which do not break chiral symmetry. 
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value of ~ 9, which would give ujq{p) ~ 3T, or ~ 1.2 GeV. Furthermore, their 
rUq are calculated with bare quark masses which increases rriq, whereas our 
calculations are in the chiral limit. At least for the moment, our use of = 1 
GeV seems not unreasonable. 

In spite of such large uncertainties in the mass, we will be able to carry out 
the calculations of meson binding, expressing it in units of mq{T), so that 
its precise value won't matter. In particularly, we will show that at T — > 
Tc{unquenched) the binding of the pion and sigma mesons will go to — 2mq(T), 
making them massless. The spin-spin and instanton molecule interaction will 
at this point be less attractive for vector and axial mesons, which will thus 
have a finite limit ^ of Mv^A/mq at T Tc. 

1.3 The "hard glue" and instanton molecules atT>Tc 

As we already mentioned in the introduction, the trace of the stress tensor 
T^ij, — € — 3p, is not "melted away" at T > Tc, but is instead only a factor 
2 smaller than at T = 0. Nice details on this observable have been recently 
provided by David Miller [23], who used a set of the Bielefeld group lattice 
data and showed that the hard glue is only melted by T ^ 400 MeV. 

The natural explanation of the existence of such "hard glue" or "epoxy" com- 
ponent, which survives chiral restoration and dcconfinement, was given by the 
theory of instantons. It has been argued [16] that chiral restoration is not re- 
lated to the instanton suppression at T > Tc, as was previously thought, but to 
their rearrangement, from the quasi-random ensemble to that made of corre- 
lated instanton-anti-instanton pairs (or other clusters with net zero topological 
charge). Some details of the quantum/statistical mechanics of a // molecules 
can be found in [24], the finite T simulations of instanton ensemble has been 
shown in [25], for other references see [15]. 

For a single molecule the contribution to the partition function can be written 
as 



where fl, fl are 12-dimensional collective variables for 1,1, namely size, 4 po- 

^ Note however that the Harada-Yamawaki vector manifestation (VM) [21] of chiral 
symmetry predicts that the vector mesons also reach zero mass at from below 
in accordance with the BR scaling. In a forthcoming publication [22], an argument 
wiU be developed to the effect that the vector mesons also go massless as T — > 
from above. 




(2) 
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sitions and 7 SU{3)c angles^. Here Tjj is the matrix element of the Dirac 
operator between 2 zero modes, that of the instanton and anti-instanton. 

Instead of evaluating the coupling constant from first principles, one can do in 
phenomenologically, in a NJL-like framework. Brown et al [17] (BGLR) showed 
that the scaled attraction which breaks chiral symmetry below Tc gives way to 
one only slightly (~ 6%) lower in which chiral symmetry is restored at Tc (un- 
quenched). Whereas the Nambu Jona-Lasinio (NJL) below Tc was connected 
to the "soft glue" part of (which may be associated with the spontaneously 
generated part of ^ ) the scale anomaly which gets restored with chiral symme- 
try, at least half, if not more, of the glue - which we called "epoxy"- remained 
at Tc and for some distance above, being melted only gradually with increasing 
T. The "hardness" of this glue explains why Tc{quenched) is so much higher 
(by nearly 50 %) than Tc{unquenched) . 

1.4 The "screening" masses and states 

The issue was raised first by DeTar et al [18] and explained by Koch et al.[30]. 
The "screening" states are formed due to a specific nonperturbative phenom- 
ena in the magnetic sector of high-T QCD related with "dynamical confine- 
ment of spatial Wilson loops". These phenomena exist at all T > Tc and 
are best explained [30] in a "funny space" in which the coordinates t and z 
had been interchanged. In this way the new "temperature" in the old z di- 
rection was zero, whereas the new z is compactified to (ttT)"^, essentially a 

^ One of the angles conjugated to the Gell-Mann matrix As, does not rotate the 
instanton solution since it has only 2 colors. 

^ An elaboration is perhaps in order on the nature of symmetry breaking involved 
here with the conformal invariance. As mentioned above and elsewhere [17,26], the 
trace anomaly of QCD associated with the scale symmetry breaking has, roughly 
speaking, two components: One "soft" component which is locked to chiral sym- 
metry and the other "hard" or "epoxy" component which is not directly tied to 
chiral symmetry. In fact in [26], it was explicitly shown how in dense medium the 
soft component "locks onto" the property of chiral symmetry, with the melting of 
the soft component corresponding to the melting of the quark condensate. This is 
the notion used in the early discussion of BR scaling [27]. The simplest way to 
understand this phenomenon in the case at hand is to consider the soft component 
as resulting from an "induced (or spontaneous) symmetry breaking" and the hard 
component as an explicit symmetry breaking of conformal symmetry. As is known 
since a long time [28], conformal symmetry can be spontaneously broken only if 
it is also explicitly broken. Thus we can identify the soft component of the glue 
that melts across the phase transition as arising from spontaneous breaking in the 
presence o/of an explicit breaking which remains intact across the phase boundary. 
This is analogous to "induced symmetry breaking (ISB)" of Lorentz symmetry in 
the presence of chemical potential discussed in [29] . 
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dimensional reduction. 



In the Koch et ah work[30], in the old time (new z direction), the Coulomb 
potential came from a periodic array of charges, because of the periodicity in 
this direction, which in the large-T" limit became that from a wire 

= ln(rx/2/?) (3) 

where P — T~^. Such a potential had earher been obtained in QCD by Hansson 
& Zahed [31]. 

A linearly rising potential from the space-like string tension was added, be- 
cause it was felt that confinement on top of the logarithmic attraction was 
needed to hold the quark and anti-quark together. The —g'^{di ■d2/r) interac- 
tions were, however, in the two-dimensional x, y directions in the disc through 
current loops, and give about one eighth of the attraction we shall find. An- 
other eighth, as we outline below, comes from the usual Coulomb attraction 
which goes as —g^/r at larger distances, and which is sufficient to bind the 
quark and anti-quark. The largest part of the binding around T^. comes from 
the instanton molecule interaction. 

Koch et al.[30] ignored the effective Coulomb interaction (in the 2d form of 
that from a charged wire) because of the dominance of the space-like string 
tension. Putting together the above forces produced the two-dimensional tt 
and p wave functions in the x, y-planc. Since the x and y directions were 
unchanged, on this plane wc should recover these same wave functions in 
a projected four- dimensional QCD calculation, but with increased binding 
because of the additional interactions included. 

The important feature of DeTar's spatially propagated states was that quarks 
were still confined in colorless states above (unquenched) [18]. However in 
our work we discuss real propagating states, and so in the deconfined phase, 
T > Tc, those can be colored. The non-singlet qq states are only the color octet 
ones, which is a channel with color repulsive force and is obviously unbound: 
but states made of qg and gg type have colored channels with an attraction 
and should exist. One more famous example of colored bound states is the 
qq Cooper pairs, leading to color superconductivity phases [32]. Although we 
could have discussed all of them in the same way as we did the qq ones, we 
defer this discussion to future works. The 32 lowest qq states we consider 
here are, however, colorless. ^ Thus in the temperature region up to ~ 2Tc 

^ In total, there are 32 degrees of freedom of the qq bound states which He in 
mass well below the masses of quarks and gluons, and are therefore the relevant 
variables for the thermodynamics. This is somewhat fewer than the 37 equivalent 
boson degrees of freedom arising from perturbative quarks and gluons. This may be 
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(unquenched), essentially up to the temperatures reached by RHIC, we have 
a dynamical color confinement. 



2 Binding of the qq states 

2.1 The Coulomb interaction and the relativistic effects 

At T > Tc the charge is screened rather than confined [20] , and so the potential 
has a general Debye form 



(Note that we use a (somewhat nonstandard) definition in which as absorbs 
the 4/3 color factor.) The general tenet of QCD tells us that the strength of 
the color Coulomb should run. We know that perturbatively it should run as 



with Aqcd ~ 0.25 GeV. The issue is what happens when the coupling is no 
longer small. In vacuum we know that the electric field is ultimately confined 
to a string, producing a linear potential. 

In the plasma phase this does not happen, and SZ assumed that the charge 
runs to larger values, which may explain the weak binding at rather high T we 
discussed in the introduction. Lattice results produce potentials which, when 
fitted in the form V{r) — — Aexp(— mr) + B with constant A,B indeed in- 
dicate ^ that A{T) grows above untill its maximum at T = 1.4 before 
starting to decrease logarithmically at high T. The maximal value of the aver- 
age coupling max{A) ^ 1/2. This is the value which will keep charmonium 
bound, as found by Asakawa and Hatsuda, up to 1.6Tc[6]. 



a partial explanation of the ~ 20% lower number of degrees of freedom than would 
result from quarks and gluons, found in lattice calculations. 

^ We thank Stefano Fortunato for useful discussion of the issue and for the fits he 

provided to us. 

The maximal value of the running coupling in SZ was taken to be 1, but its value 
at the most relevant distances is about 1/2 also. 
We thank Loic Grandchamps for this calculation. 




(4) 
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(5) 



log(g/AQCD) 
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Running of the coupling is not very important for this work in which we are 
mostly interested in deeply bound states related with short enough distances. 
Therefore we will simply keep it as a non-running constant, selecting an ap- 
propriate average value. 

It is well known in the point charge Coulomb problem (QED) that when Za 
is increased and the total energy reaches zero there is a singularity, preventing 
solutions for larger Za. In the problem of the "sparking of the vacuum" in 
relativistic heavy ion collisions, the solution of the problem was found by 
approximating the nuclei by a uniformly charged sphere; for a review of the 
history see Rafelski et al. [33]. As a result of such regularization, the bound 
electron level continues past zero to — m, at which point e^e^ production 
becomes possible around the critical value of Z^.^ = 169. In short, the problem 
of the point Coulomb charge could be taken care of by choosing a distributed 
electric field which began from zero at the origin. 

In QCD the charge at the origin is switched off by asymptotic freedom, the 
coupling which runs to zero value at the origin. A cloud of virtual fields making 
the charge is thus "empty inside" . We will model a resulting potential for the 
color Coulomb interaction by simply setting the electric field equal to zero at 
r = 0, letting it decrease (increase in attraction) going outward . We can 
most simply do this by choosing a charge distribution which is constant out 
to R, the radius of the meson. If the original 2mq mass were to be lowered 
to zero by the color Coulomb interaction and instanton molecule interaction, 
then the radius of the final molecule will be 



although the rms radius will be substantially greater with the instanton molecule 
interactions playing the main role around T^. 



V^-as^is-^], r<R 
2R\ R^j' 

= r>R. (7) 

r 

This V has the correct general characteristics. As noted above, the electric field 
E must be zero at r = 0. It is also easy to see that V must drop off as r^/ R? as 



Just at r = the quark and anti-quark are on top of each other, so the electric 
field is clearly zero. (Although the quark and anti-quark are point particles, their 
wave functions will be distributed.) For two rigid spheres, V would take up the 1/r 
behavior only at 2i?, but there will be some flattening. It will become clear that our 
main conclusions are independent of these details. 



10 



the two spheres corresponding to the qiiark and anti-quark wave functions are 
pulled apart. Precisely where the potential begins the 1/r behavior may well 
depend upon polarization effects of the charge, the + and — charges attracting 
each other, but it will be somewhere between R and 2i?, since the undisturbed 
wave functions of quark and anti-quark cease to overlap here. 

The qq system is similar to positronium in the equality of masses of the two 
constituents. Since the main term value is^^ the 4, rather than 2 in 

hydrogen, coming from the reduced mass, one might think that the Coulomb, 
velocity- velocity and other interactions would have to be attractive and 8 times 
greater than this term value in order to bring the 2m.g in thermal masses to 
zero. However, this does not take into account the increase in reduced mass 
with a. Breit and Brown [34] found an a^/4 increase in the reduced mass 
with a, or 25% for ct = 1, to that order. It should be noted that in the Hund 
and Pilkuhn [35] prescription the reduced mass becomes /i — rrig/E, which 
increases as E drops. 

We first proceed to solve the Coulomb problem, noting that this gives us the 
solution to compare with the quenched lattice gauge simulations, which do 
not include quark loops. 

Having laid out our procedure, we shall proceed with approximations. First 
of all, we ignore spin effects in getting a Klein-Gordon equation. The chirally 
restored one-body equation which has now left-right mixing is given by 

{po + d-p)7/;^0. (8) 
Expressing ip in two-component wave functions $ and ^, one has 

Po* = -<? • P^: (9) 
giving the chirally restored wave function on ^ 



(po-a- p—a ■ p W = 0. (10) 
V Po J 

Here 



Po^Ev^E + as/r. (11) 



To make comparison with QCD we should use Z'^a rather than a and remember 
the Casimir factor 4/3. 
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Neglecting spin effects, a-p commutes with po, giving the Klein-Gordon equa- 
tion Pq — = ^- We now introduce the effective (thermal) mass, so that the 
equations for quark and hole can be solved simultaneously following [35]; 



[e-V{r)f - -p'' V'(r) = 



(12) 



where p is momentum operator, and the reduced energy and mass are e = 
{E"^ — m\ — m|)/2£', = mim2/ E with mi = m2 = ruq. 



Furthermore from eq.(9). 



(a) = (*^ a^) + ($^ am = ^- —i[a x p\). (13) 

Po Po 

If a is parallel to p, as in states of good helicity, the second term does not 
contribute. Prom the chirally restored Dirac equation (8), ignoring spin effects 
such as the spin-orbit interaction which is zero in S-states we are considering, 
we find p^ — p^. 

Brown [36] showed that in a stationary state the EM interaction Hamiltonian 
between fermions is 



e2 



^int = — (1 - ai • 02) , (14) 

where the di^2 are the velocities. Applying (14) to the chirally restored domain 
of QCD, we expect 



ri\r 



'int 



for oil ■ a2 = —1 
for di- Q.2 = +1 



(15) 



2.2 The spin-spin interaction 



The nonrelativistic form of the spin-spin interaction, in the delta- function 
form, may give an impression that it is maximal at the smallest distances. 
However this is not true, as becomes clear if the relativistic motion is included 
in full, and in fact at r ^ it is suppressed. At large r, when particle motion 
is slow, it is of course again suppressed, thus contributing mostly at some 
intermediate distances. 

This fact is clear already from the derivation of the Is-state hyperfine splitting 
Fermi-Breit due to hyperfine interaction in hydrogen from 1930 [37] given by 
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5H^^-ia.fl) fd^M^ ^ (16) 

Note the complete denominator, which non-relativistically is just substituted 
by m alone, but in fact contains the potential and is singular at r — 0. The 
derivative of the e^/r in the denominator insured that the electric field was 
zero at r = 0. Here a is the electron spin, /I the proton magnetic moment. In 
eq (16) the derivative can then be turned around to act on 'j/'V; ^-^id to order 
a—1 and with the e^/r neglected in the denominator, one has 



5/7^-^(a.^)-^V'^(0), (17) 

with 'ip taken to be the nonrelativistic Is wave function to lowest order in a. 

The hyperfine structure is obtained by letting the first p in eq. (10) act on the 
Po ^ and the second p go p+ -s/al A with 



with /I the magnetic moment of the anti-quark. One finds that the hyperfine 
structure is [37] 



H^(s = ^V^s^-[ExA] (19) 
Po 

— * 

where E is the color electric field. Thus, 



^ ^ _ y/al \E\ / a ■ pi ^ a -f fl-f \ _ 2 \E\a ■ jl 



where |^| = 2as/r^- As in the hydrogen atom, the magnetic moments of 
quarks and anti-quarks are 



= T (21) 

except that the Dirac mass niq^q = and po, in which the potential is increased 
by a factor of 2 to take into account the velocity- velocity interaction, is now 



Po^E + 2{as/r) (22) 



13 



for QCD so that in terms of the quark and anti-quark magnetic moment 



operators , 



2 

^^hfs = - ^ 2 (1^1 ■ Z^?-) • (23) 

o Par 

Of course, our po for the chirally restored regime has substantial r dependence 
(whereas the e/r in the hydrogen atom is generally neglected, and + m is 
taken to be 2m, so that /ig = — e/2me). From Fig. 3 it will be seen that (square 
of) the wave function is large just where cts/r is large. 

For rough estimates we use averages. We see that, as in Table 1, if £■ is to be 
brought down by ~ O.Sm^ for the a and tt by the Coulomb interaction, then 



1 1 

2(a,/r) ~ -rUg ~ -po (24) 
so that with ~ 0.5, 



{r-') ~ ^m,. (25) 

We next see that this is consistent with the spin splitting forming a fine struc- 
ture of the two groups, the lower lying a and tt, and the slightly higher lying 
vectors and axial vectors. Using our above estimates, we obtain 



so that for the a and tt where iJi • (T2 = — 3 we have 



Tfl 

{Hnfs) - (27) 

the approximate equality holding when ag — 0.5. Note that the hyperfine 
effect is negligible for the ag ~ 0.5. Although formally eq. (26) looks like 
the hyperfine structure in the chirally broken sector, it is really completely 
different in makeup. 

In our expression for (Hh/s) we have the r dependence as {pQr)~'^r~^ and 
Por — 4, basically because the Coulomb interaction lowers the tt and a only 



The electric field is denoted as E which should be distinguished from the energy 

E. 
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1/4 of the way to zero mass. This explains most of the smallness of the spin- 
dependent interaction. 

A recently renewed discussion of spin-spin and spin-orbit interactions in a 
relativistic bound states has been made Shuryak and Zahed [38], who derived 
their form for both weak and strong coupling limits. Curiously enough, the 
spin-spin term changes sign between these two limits: perhaps this is another 
reason why at intermediate coupling considered in this work it happens to be 
so small. 



2.3 The effect of the instanton molecules 



The effective interaction is calculated as follows. The propagator is written as 



S(.,!,) = E%^ (28) 
^ X-\-im 

where A, 0a are eigenvalues and eigenfunctions of the Dirac operator. If the 
gauge field configuration is the anti-instanton-instantopn (//) molecule, the 2 
lowest eigenvalues are A = ilT/^l and their eigenfunctions are simple combi- 
nations of zero modes for instanton and anti-instanton 0/ ± 0/. The sum of 
those leads to 



S{x, y) = -^(ct>\{x)ct>j{y) + ct>\{x)Mv)) (29) 

This can be interpreted as follows: 0/(a;) is the amplitude to go from point x 
to the instanton 7, the other phi is the amplitude to appear from the anti- 
instanton, and the factor in front is the propagator between I and /. 

For propagation of 2 quarks, say of opposite fiavors, one can then draw the two 
diagrams shown in Fig.l. The diagram (a) read from left to right has quarks 
of opposite chirality, so it contributes to scalar and pseudoscalar mesons; the 
diagram (b) has the same chirality of quark and anti-quark, so it contributes 
to vector and axial vector channels. One can see from the figure that in the 
former case both qq go into the same instanton, while in the latter they have 
to go to the opposite ones. As a result, the locality of the former vertex is 
given by the instanton size p and of the latter by the size of the molecule. 

If all 4 points are far from the instanton and anti-instanton, (j){x) ~ l/x^ mod- 
ulo a constant spinor. This power of the distance corresponds to free prop- 
agator for a massless quark, 5*0 = l/(27r^x^), which should be "amputated". 
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Fig. 1. The diagrams generating the effective Lagranfian 

leading to a constant vertex function, or effective Lagrangian [25]. The ef- 
fective interaction between quarks is conveniently calculated by rearranging 
the exchange terms into a direct interaction. The resulting Fierz symmetric 
Lagrangian reads [25] 



Lmolsym = G { [{i^r'^i^f - (V^T^TsV')' 



1 r 



(VJt^t^V')^ + (V^r«7^75V^)'] + (^^7^75^)'} + i^s, (30) 



with the last complicated term containing color-octet qq pairs only. The cou- 
pling constant is now proportional to the density of molecules 



G = j^j n{puP2) dpidp2 ^(27rpi)2(27rp2)' ■ (31) 

If the molecules are incompletely polarized in the (Euclidean) time direction 
then all vector terms are modified accordingly, because {'ijjjijrip) ~ Z/^, the 
only vector available. 

We start by explaining the issue of the "polarization" of the // molecules. 
Let us put / at the origin and in the standard SU(3) orientation without 
rotation. Let us introduce the 4-d polar angle 64: the position vector of /, 
called Zfj, is such that Z4 = cos{94)^J^. We will refer to the 64 = case as 
molecules completely polarized in the time direction: as discussed in [25] this 
is where the maximum of the Z is. The very robust maximum corresponds to 
"half Matsubara box" , Z4 = 1/ (2T), and if 64 is nonzero Z is smaller, basically 
because of large distances in the T/j. Let us estimate the effect of that, ignoring 
the gauge action. 



which can be expanded into Gaussian form at small 64. For 2 flavors and 
T — Tc ^ l/(4p), the root mean square polarization angle is (6*4) ~ 0.55. We 
will need below 
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cos^aei)) ^ .72 



siri^ae^)) ^ .28 



(33) 



corresponding to (6*4) = 32°. 

Let us remember that for vector particles propagating in matter one can de- 
fine in general 2 structures in the polarization operator, the longitudinal and 
transverse ones IIl,t, which are related to 1 longitudinal and 2 transverse 
modes. Those are related to Cartesian components as follows 



Hoo = Ul, Hon = ^ (34) 

TT fx PmPn\ -f-f I PmPn -pr /qc:\ 
= l^mn ^ I H ^ — 11^, (35) 

where the Latin indices are 1-3 and Greek 0-3, cu — pq. This polarization tensor 
satisfies the conservation law P/^U/^i,, eliminating the 4-th component. 

For complete polarization of molecules, 64 = and only the zeroth component 
(the longitudinal Ul component) is coupled, while in general the coefficients 
include sin^ and cos^ of the angles determined above (33) in 11^ and 11^,, 
respectively. 

This results in additional nonlocal correction factors, which for the scalar- 
pseudoscalar channels is 



^S,PS 
nonlocal 



x(o) 



(36) 



and for the vector-axial-vector ones 



(37) 

where the additional angular bracket in the latter case comes about from 
averaging over the molecular orientation relative to the time axes. Note that 
in the former well as in the latter for ^4 = 0, the correction factors 

are 1, due to the normalization condition / d'^x\(j)i{x)\'^ = 1, for weakly bound 
(large size) states, for which the factor x(a;)/,x(0) can be approximated by 
1 and taken out of the integral. We treat the ratios (x(a;)/x(0))^ = F as a 
vertex correction in our calculations in the next section, and estimate their 
size in the Appendix. 

The effective interaction of light quarks, the Fierz symmetric instanton molecule 
Lagrangian of Schafer et al. [25] in Minkowski space, is of the form 



nonlocal 
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Fig. 2. Looking down on the instanton molecule close-packed around the antiperiodic 
time direction at Tc- Here p is the | fm radius of the instanton or anti-instanton 
and r = 1/4/3 ~ 150 MeV. Just at Tc the instanton molecule is fully polarized in 
the r direction. The lines going in and out of the instanton and anti-instanton are 
the quark zero modes. The outer quark zero modes are broken open to give the 
quark-quark interaction which forms the instanton induced interaction. 



+ 4 



^ (^r"7^^)2 - (^r"7M75^)' (38) 



The value of G at Tc was obtained by BGLR [17] as 



G = 3.83 GeV-^ (39) 

Here = (r, 1) is a four-component vector. This Lagrangian gives attraction 
in the a, 5, tt, p and Ai sectors. We will find that the mesons are bound by both 
the color Coulomb and the four-Fermi instanton molecule interactions. 



We should explain how the chirally broken NJL is relayed into the chirally 
restored NJL as T goes upwards through T^. The 't Hooft instanton-driven 
interaction has been included in the chirally broken NJL, and undoubtedly 
enters into the interactions which bring the m^r to zero and the to 2mg, 
although their role relative to the other interactions is not clear in this do- 
main. As noted in BGLR[17], about half of the total bag constant (confor- 
mal anomaly) coming from the spontaneous breaking of chiral symmetry (the 
restoration of which is associated with Brown/Rho scaling) is transferred from 
that of the binding energy of the negative energy nucleons, which are the rel- 
evant T = 0, p = variables to the rearrangement of the random instanton 
vacuum into the instanton molecules (which do not break chiral symmetry) 
and the other half goes into melting the soft glue as the nucleons loosen into 
constituent quarks. 

BGLR[17] showed that NJL could be continued from the chirally broken region 

As often the case in condensed matter and particle physics, attractive four-Fermi 
interactions could - via anomalous dimensions - figure crucially in phase transitions, 
so may play a more important role near T^- This may be particularly relevant to 
other mesons than the tt and a. It is shown in [22] with a schematic model that the 
vector mesons do undergo a phase change very similar to that of the tt and a. 
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up to the chirally restored one with the ~ 6% decrease in the NJL couphng 



(7 = 4.08 GeV 
(7 = 3.83 GeV 



-2 



chirally broken phase 
chirally restored phase. 



-2 



(40) 



We suggest that the fragility of chiral symmetry breaking as T moves upwards 
to Tc signifies the importance of the instanton molecules in this region of 
T. With only an ~ 6% decrease in G, the chiral restoration transition is 
effectuated. 

Whereas in the chirally broken NJL, the vector interactions are about the same 
size as the scalar ones, one can see from the instanton molecule Lagrangian, 
eq. (38), that the p and Ai interactions are a factor of 4 smaller than the tt 
and a ones. However, as we have outUned, taking the instanton molecule to be 
polarized along the time axis, the time components of p and Ai interactions 
are built up a factor of 4, at the expense of the spatial ones. Thus, in the 
classical approximation, neglecting fluctuations in ^4, the vector and axial 
vector modes would be degenerate with a and tt. 

For the local 4-fermion interaction with the coupling constant G the energy 
shift is given simply by 



More generally, for the non-local interaction induced by molecules, one should 
project on the (2-body) wave function of the bound state ^(a;, y)) = exp[iPn{x+ 
y)^/2]x{x — y) where, for the stationary state, there is no dependence on rel- 
ative time xq — yo. 

The interaction is attractive in all channels, so that if the Coulomb interaction 
does not bring the meson mass all the way to zero, the instanton molecule 
interactions will, in particular near the phase transition temperature Tc, where 
the wave functions are very compact. 

In the NJL-like instanton molecule Lagrangian eq. (38) the interactions are 
four-point in nature. We can convert these to an instanton molecule interaction 
by constructing a pseudo-potential V = C6{x) for the ipip — > ijjilj scattering 
amplitude. Here {■ijj5{x)ijj) = ^^(0). Therefore G — Gtjjijj{0) with G the inter- 
action of eq. (39), gives the strength of the pseudo-potential. 

Let us mention some estimates for the mass of the vector/axial mesons. Keep- 
ing as above the thermal quark mass at niq = IGeV and the as = 0.5 results 
from Table 1, we will see that the Coulomb binding of ~ 0.5 GeV plus that 



SE = G\iP{0)\^ 



(41) 
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due to instanton molecules ~ 1 GeV is able to make tt, a mesons massless, 
once we sum loops. 

In the Appendix we estimated that the corrections for nonlocality were roughly 
equal for tt, a and p, Ai, so that the masses of the latter scale as cos^((^4)) 
and sin^((^^4)) relative to the 2mq ~ 2 GeV binding energy of the tt and a. 
From our estimates Eq. (33) we see that the (predominantly) longitudinal p, 
Ai masses will come at ~ 560 MeV and the transverse p, A masses at 1440 
MeV; i.e., at roughly the free p and Ai masses, respectively. However, the 
quasiparticle composed of coupled p and 2ti components ( "rhosobar" ) may lie 
lower in energy when the interaction is diagonalized. 

2.4 The resulting qq binding 

We first construct the bound states for T ^ Tc{unquenched) , at temperature 
close enough to Tc so that we can take the running coupling constants at 
T = Tc + e. The fact that we are above Tc is important, because the A^sb ~ 
4vr/^ ~ 1 GeV which characterizes the broken symmetry state below Tc no 
longer sets the scale. Until we discover the relevant variables above Tc we 
are unable to find the scale that sets as = ||^, the color Coulomb coupling 
constant. 

Following SZ [11], we adopt quark- ant i- quark bound states to give the rel- 
evant unperturbed representation and, the instanton molecule gas [25] as a 
convenient framework. In particular, Adami et al. [39], Koch and Brown [40], 
and BGLR [17] have shown that ^ 50% of the gluon condensate is not melted 
at T = Tc. The assumption motivated by Ilgenfritz & Shuryak [16] is then 
that the glue that is left rearranges itself into gluon molecules around T = Tc, 
i.e., what BGLR call "epoxy". We have quantitatively determined couplings 
for the mesons in the instanton molecule gas by extending the lower energy 
NJL in the chiral symmetry breaking region up through Tc [17]. We set these 
couplings in order to fit Miller's [41] lattice gauge results for the melting of 
the soft glue. 

In Fig. 3 we show that if we choose — 1 (effectively c^s = 2 by the doubhng 
in Eq. (15)) as would be required to enter the strong coupling region considered 
by Shuryak and Zahed[14] we bring the meson mass down by —1.36 rriq from 
their unperturbed 2 niq. However, we switch to the region of ~ 0.5, which 
is required by charmonium (intermediate coupling). In Table 1 we summarize 
the Coulomb binding for a few choices of as- 

In the case of the instanton molecule interaction the coupling constant G = 
3.83 GeV^^ is dimensionful, so that its contribution to the molecule energy 
scales as G m^. (Since we take as — 0.5 and will find that with inclusion of 
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R 




Fig. 3. The color Coulomb potential V and the corresponding wave function ip for 
relativistic Klein-Gordon case. The interaction eq. (7) with R = h/2mq was used. 
The ground state energy with ag = I corresponds to -Egj-ound = 0.645 niq. The 
minimum of the potential at the origin is at —3mq here. 

the velocity- velocity interaction the effective ctg will be 1, powers of a will not 
affect our answer. We will use ruq — 1 GeV, essentially the lattice result for |Tc 
and 3Tc[13], which works well in our schematic model.) Of course, in QCD the 
Polyakov line goes to zero at Tc, indicating an infinite quark mass below Tc, 
i.e., confinement. Just at Tc the logarithmically increasing confinement force 
will not play much of a role because the dynamic confinement holds the meson 
size to ~ h/mqC, or ~ 0.2 fm with our assumption of niq = 1 GeV. (Later we 
shall see that the rms radius is ~ 0.3 fm.) Since we normalize the instanton 
molecule force, extrapolating it through Tc{unquenched) , and obtain the color 
Coulomb force from charmonium, our niq is pretty well determined. However, 
our niq — 1 GeV is for the unquenched system and at a temperature where 
the instanton molecules play an important role. 

Given these caveats, we may still try to compare our Coulomb result with the 
lowest peak of Asakawa et al.[42] which is at ~ 2 GeV for T = lAT^ ~ 0.38 
GeV and for Petreczky at ^ 5T ~ 2.030 GeV for T = 1.5T^ ~ 0.406 GeV 
where wc used the Asakawa et al. Tc (quenched). We wish to note that: (i) 
These temperatures are in the region of temperatures estimated to be reached 
at RHIC, just following the color glass phase (which is estimated to last ~ 1/3 
fm/c). Indeed, Kolb et al. begin hydrodynamics at T = 360 MeV. (ii) These 
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Ots 


A£^Couiomb [GeV] 




A-B4-pomt[GeV] 


0.50 


—0.48 


0.36 


—0.99 


0.55 


-0.60 


0.31 


-1.39 


0.60 


-0.71 


0.28 


-1.83 


1.00 


-1.36 


0.14 


-7.57 



Table 1 

Binding energies from color Coulomb interaction and the corresponding rms radii 
for various (effectively, 2as including velocity- velocity interaction). 4-point in- 
teractions are calculated using the parameters obtained from color Coulomb inter- 
action. 




Fig. 4. Coulomb molecule. The wavy line on the left represents the momentum 
transfer necessary to produce the molecule. The double line denotes the Furry rep- 
resentation, i.e., the full propagator in the Coulomb potential. 

are in the region of temperatures estimated by SZ[11] to be those for which 
bound mesons form. 

We are unable to extend our consideration to higher temperatures, where the 
situation may move towards the perturbative one, but we believe that lattice 
calculations do support our scenario that the QGP contains large component 
of bound mesons from T ~ 170 MeV up to T ~ 400 MeV. 

For Us = 0.5, which is the value required to bind charmonium up through 
T — 1.6Tc, we find that the Coulomb interaction binds the molecule by ~ 0.5 
GeV, the instanton molecule interaction by ~ 1.5 GeV. However, the finite 
size of the ip^ip of the instanton zero mode could cut the latter down by an 
estimated ~ 50% (See Appendix). As in the usual NJL, there will be higher 
order bubbles, which couple the Coulomb and instanton molecule effects. 
We draw the Coulomb molecule in Fig. 4, where the double lines denote the 
Furry representation (Coulomb eigenfunction for quark and anti-quark in the 
molecule) . 

The four-point instanton molecule interaction is shown in Fig. 5. There will 
be higher-order effects as shown in Fig. 6, of the 4-point interaction used 
in higher-order between Coulomb eigenstates which always end in a 4-point 
interaction. The energy of the propagators has been lowered from the 2 GeV 

Alternatively, the (^-function of the 4-point interaction should be included in the 
Klein-Gordon equation, where it will change the energy. 



22 



Fig. 5. The four-point instanton molecule interaction between Coulomb eigenstates. 
The {'ipipy intersect at the thick point. 




Fig. 6. Higher order effects of four-point interaction. 

of the two noninteracting quarks to 1.5 GeV by the Coulomb interaction. The 
series beginning with terms in Figs. 4— 6 is 

^ -0.5GeV - IGeV ^ - - + . . . 

l.SGeV (1.5GeV)2 

= -0.5GeV--^^J^. (42) 

1.5GeV 

Now AE — —1.25 GeV is accomphshed for F — 0.5, a reasonable assumption 
from the estimate in the Appendix. 

Working in the Furry representation (including the Coulomb potential ex- 
actly), we have a —0.5 GeV shift already from the Coulomb wave functions. 
This means that we must obtain AE — —1.5 GeV to compensate for the 
2mq — 2 GeV, in order to bring the tt and a masses to zero. The four-point 
interaction is a constant, at a given temperature, so this problem is just the 
extended schematic model of nuclear vibrations (See Sec. V of Brown [44], 
where simple analytical solutions are given). 

Our eq. (42) corresponds to the Tamm-Dancoff solution, summing loops going 
only forward in time. If AE decreases —0.75 GeV in this approximation, then 
when backward going graphs are added, then AE will decrease by twice this 
amount [44], or the —1.5 GeV necessary to bring the tt and a energy to zero. 
(There is an analogy between the treatment of the spurion and the Goldstone 
mode in hadronic physics.) Of course, forward and backward going loops are 
summed in the Bethe-Salpeter equation to give the NJL in the broken sym- 
metry sector, but the actual summation is more complicated there, because 
the intermediate state energies are not degenerate. 

In detail, with our estimated F = (0.75)^ from the Appendix and the 4-point 
energies from Table 1, our n and a excitations without inclusion of backwards 
going graphs are brought down 58% of the way from —0.5 GeV to —2 GeV; 
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i.e., slightly too far. We have not made the adjustment down to 50%, because 
the uncertainties in our estimate of F, etc., do not warrant greater accuracy. 



2.5 Comparison with Lattice Calculations 

Lattice calculations propagating quarks and anti-quarks in the spatial direc- 
tion gave evidence of DcTar's dynamical confinement [45]. The dynamical 
confinement comes about from the attractive agC^i ■ 012/ or current-current 
interaction in our eq. (14). This is the same in the "funny space", since the x 
and y directions, in which the current loops lie, are as in the real space. Thus 
we find from eq. (14) that the introduction of the Coulomb interaction dou- 
bles the effects from the current-current interaction. The instanton molecule 
interaction gives a factor of several times the current- current interaction. It 
should be included in the unquenched lattice calculations. 

It is instructive to examine the results of Bernard et al. [45] (reproduced in 
Fig. 8. 6 of Adami and Brown [46]) for two-dimensional wave functions at 
T = 210 and T = 350 MeV of the chirally restored tt and p. These wave func- 
tions are given in physical units. They have been measured on the (x, ?/)-plane 
and we think of them as projections onto this plane of the three-dimensional 
wave function. The pion wave function at T = 210 MeV is seen to drop off ex- 
ponentially, decreasing to 0.02 of its r = value by r = 1.15 fm. With a wave 
function C exp(— /tr), this indicates a k of 3.4 fm~^ = 0.68 GeV. This gives an 
rms radius for the two-dimensional wave function of v^/2/t=1.8 GcV-^ = 0.36 
fm. For the three-dimensional one we would multiply by \JS/2 ~ 1.22; thus 
one arrives at an rms radius of ~ 0.44 fm. This is not much larger than the 
0.36 fm for our chosen a — 0.5, and it should not be because of the relatively 
small role played by the Coulomb term. 

It can be seen that the p-meson wave function does not drop exponentially 
until 2 — 3 fm, showing it to be larger in extent than the pion. It then appears 
to decrease somewhat more slowly than the pion wave function, although the 
errors are such that there is uncertainty in this. The true quasiparticle above 
Tc may, however, be a linear combination of p and 2 pion states. 

The instanton molecule model, as the DeTar dynamical confinement, involves 
an analytical continuation from imaginary to real time, although it is quite 
simple in the former. The SZ [11] Coulomb-bound gas of mesons is formulated 
directly in real time. It is simple and straightforward. We view it to be useful 
to give our alternative instanton molecule formulation, because it enables us 
to make contact with the lattice calculation, especially with those of the gluon 
condensate. 
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Hadronic spectral functions above the QCD phase transition have been calcu- 
lated in quenched lattice gauge simulations by Asakawa, Hatsuda, & Nakahara 
[42] and by Petreczky [43]. 

In the quenched approximation, the wavy line in Fig. 4 can be interpreted as 
the source current which creates a valence quark and anti-quark. These are 
allowed to exchange any number of gluons, so it is appropriate to let them 
exchange any number of Coulomb interactions, as well as develop thermal 
masses; i.e., the Coulomb problem is that which we outlined. 

We differ in detail from SZ[11], in that we consider the 32 normal modes of qq 
states, in which the quark and anti-quark, or quark and hole, have opposite 
helicities so as to benefit from the current- current interaction. SZ have the 
bound states of quark and anti-quark, or quark and hole, and, separately, 
those of gluons, the glueballs. Because of the larger Casimir operator for the 
gluons, these bind at a higher temperature than the quarks, so there would 
be two regions of temperature in which the molecules break up. The number 
of degrees of freedom, 40 quarks and 16 gluons, once they break up, would be 
different. Presumably they would have the same qq bound states for T 
as we, but, in addition, would have the glueballs. 

As noted by SZ, in the region where the molecules break up, the quark ve- 
locities arc small, going to zero at zero binding energy. This means that the 
velocity-velocity interaction will be unimportant and the additional 32 de- 
grees of freedom, disfavored at lower temperatures by the velocity-velocity 
interaction, will become important, doubling the degrees of freedom. 

However, the quenched approximation would not include the quark and anti- 
quark loops of the instanton molecule interaction. Thus, at T = Tc we would 
have molecules of energy 1.5 GeV, since our Coulomb interaction gives a 0.5 
MeV binding on thermal quark and anti-quark masses, each of 1 GcV. (In fact, 
these masses were measured at 1 .5Tc{quenched) , and need not be the same at 
Tc. Furthermore, we need to mention that our Tc is that for the unquenched 
calculations, Tg = 170 — 175 MeV, because we necessarily have a situation 
with quarks and anti-quarks, especially instanton zero modes.) 

In Fig. 7 we show Asakawa & Hatsuda[6] spectral functions for (quenched) 
temperatures of 1.4Tc ~ 380 MeV and 1.9Tc ~ 515 MeV. The lower temper- 
ature is essentially that reached at RHIC following the color glass phase and 
the second temperature is higher, but probably not in the perturbative regime 
because the lattice calculations by Petreczky et al.[13] give the quark mass at 
3Tc as roughly the same as at |Tc, not increasing as gT. 

The narrow 2 GeV peak supports our identification of this temperature (T ~ 
380 MeV) being where the molecules break up, because the particle velocities 
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Fig. 7. Spectral functions of Asakawa et al.[42]. Left panel: for Nj- = 54 (T ~ 1.4Tc). 
Right panel: for A^^ = 40 (T ~ l.DTc). 

will begin from zero at breakup. In this case, the 2 GeV reproduces the sum 
of quark and anti-quark masses. We note the meson degeneracy is consistent 
with our negligible effects from spin and chiral restoration of the tt and a and 
vector and axial vectors. 

The 1.9Tc ~ 515 MeV data which should be well above the temperature at 
which the molecules break up does show definite thermal effects. As noted 
earlier, the lattice calculations give approximately the same ttIq at at 
(3/2)Tc, so we can take ruq to be constant in this region of energies. Each 
of the ~ 1 GeV quark and hole in quark and anti-quark will have a thermal 
energy of ~ 1.1 GeV, midway between the |T nonrelativistic thermal energy 
and the 3T relativistic one. Therefore, it costs ~ 4.2 GeV to produce a quark 
particle and hole or quark and anti-quark at this temperature. The half width 
of the peak is roughly consistent with the size of the proposed thermal energy. 

At low uj we probably see Landau damping in the 1.9Tc results, although at 
1.4Tc there is no such sign. This further supports our belief that the quark 
and anti-quarks are essentially at rest there. 

By rescaling the lattice calculation Peter Petreczky (private communication) 
has shown that the peaks higher than the first one are lattice artifact. 

We see that in the SZ[11] model, one can understand not only in general terms, 
but also many of the details of the Asakawa and Hatsuda lattice calculations. 



' In the LGS the zero-momentum state of the quark-anti-quark pair is projected 
out. 
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3 Conclusions 



Shuryak and Zahed have discussed the formation of the mesonic bound states 
at higher temperatures, well above Tc. They pointed out that in the formation 
of the bound state, or any one of the molecular excited states, the quark- 
quark scattering length becomes infinite, similarly for the more strongly bound 
gluon-gluon states. In this way the nearly instantaneous equilibration found by 
RHIC can be explained. As we explained in the last section, lattice calculations 
seem to support the scenario of nearly bound scalar, pseudoscalar, vector and 
axial- vector excitations at ~ 2Tc{unquenched) (~ 1.5 times the quenched Tc). 



In this work we are able to construct a smooth transition from the chirally 
broken to the chirally restored sector in terms of continuity in the masses 
of the cr and tt mesons, vanishing at T — > T^. In doing so we had to include 
relativistic effects. One of them - the velocity- velocity term related to Ampere 
law for the interacting currents - nearly doubles the effective couphng. The 
spin-spin term happens to be very small. The crucial part of strong binding 
in our picture of qq mesons (or molecules) is the quasi-local interaction due 
to instanton molecules (the "hard glue"). We found that the tight binding of 
these mesons near T^, enhances the wave function at the origin, and gives us 
additional understanding of the nonperturbative hard glue (epoxy) which is 
preserved at T > Tc- 



Thus, we believe that the material formed in RHIC was at a temperature 
where most of it is made of chirally restored mesons. Certainly this is not the 
weakly coupled quark-gluon plasma expected at high T. 



Finally, in this paper we have focused on quantum mechanical binding effects 
in the vicinity of the critical temperature coming down from above. Nice 
continuity in the spectra of the light-quark hadrons - e.g., the pions and the a 
- across the phase boundary should also hold for other excitations such as the 
vector mesons p, a;, Ai which lie slightly above tt and a because of quantum 
corrections. Since going below Tc from above involves a symmetry change from 
Wigner-Weyl to Nambu-Goldstone, there is a phase transition and to address 
this issue, it would be necessary to treat the four-fermi interactions more 
carefully than in the pseudo-potential approximation adopted here. As briefly 
noted, the true quasiparticles in the p-channel in the many-body medium may 
be linear combinations of the p found here and 27r states. 
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A Appendix: Finite Size Corrections 



The instanton is ~ 1/3 fm in radius, and the rms radius of the molecule is 
about the same for q;^ = 0.5, so there will be a correction for the finite size of 
the instanton zero modes. Their wave functions are 



ZM 



+ p2)3/2 

with p the instanton radius. Now 



(A.l) 



(x2+p2)3- 

Integrating over time 



(A.2) 



dr- 



1 



1 



^2 + (^2 + ^2)5/2 

5 



oc exp 



! exp 



2p2 




(A.3) 



Thus, we see that ■^■^ is sharply peaked, mostly lying within a radius r ~ 

^2/5 p or a volume of (^2/5 pj , or ^ 25% of the instanton molecule. This 
nonlocality will spread the initially forward peaked qq wave function over a 
volume of about 1/4 of the instanton, so we estimate F — (0.75)^ in Sec. 2.4. 

Our above estimate holds for the effect of nonlocality in the scalar and pseu- 
doscalar mesons, where its ratio to the size of the instanton is calculated. In 

the vector and axial-vector the ratio should be to that of the entire molecule, 
but the nonlocality is also over the whole molecule, and the wave function, 
due to less binding, will not be so forward peaked as in the tt and a. Thus, we 
use the same F — (0.75)^. 
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The possible considerable error in this estimate is important in determining 
the role of the vectors and axial-vectors in the thermodynamics of the system, 
but not in our main purpose of constructing continuity in the vr and a masses 
across Tc. 

Ultimately the vector and axial-vector masses may be quantitatively deter- 
mined for T ^ Tc if the lattice gauge simulation of the Asakawa and Hatsuda 
type (see Sec. 2.5) are extended to unquenched ones. 
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